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'k^ ' Abstract 

j^ ■ We investigate the infra-red stable fixed points of the Yukawa couplings in the 

minimal version of the supersymmetric standard model with R-parity violation. 
Retaining only the R-parity violating couplings of higher generations, we analyti- 
cally study the solutions of the renormalization group equations of these couplings 
together with the top- and b-quark Yukawa couplings. We show that only the B- 
violating coupling A233 approaches a non-trivial infra-red stable fixed point, whereas 
all other non-trivial fixed point solutions are either unphysical or unstable in the 
infra-red region. However, this fixed point solution predicts a top-quark Yukawa 
coupling which is incompatible with the top quark mass for any value of tan /3. 
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There is considerable interest in the study of infra-red (IR) stable fixed points of the 
standard model (SM) and its extensions, especially those of the minimal supersymmetric 
standard model (MSSM). This interest follows from the fact that in the SM (and in 
the MSSM) there are large number of unknown dimensionless Yukawa couplings, as a 
consequence of which the fermion masses cannot be predicted. One may attempt to 
relate the Yukawa couplings to the gauge couplings via the Pendleton-Ross infra-red 
stable fixed point (IRSFP) for the top-quark Yukawa coupling [|l|, or via the quasi- fixed 
point behaviour [^ . The predictive power of the SM and its supersymmetric extensions 
may, thus, be enhanced if the renormalization group (RG) running of the parameters is 
dominated by IRSFPs. Typically, these fixed points are for ratios like Yukawa couphng 
to the gauge coupling, or, in the context of supersymmetric models, the supersymmetry 
breaking tri-linear A-parameter to the gaugino mass, etc. These ratios do not always 
attain their fixed points values at the weak scale, the range between the GUT (or Planck) 
scale and the weak scale being too small for the ratios to closely approach the fixed point. 
Nevertheless, the couplings may be determined by quasi-fixed point behaviour 0, where 
the value of the Yukawa coupling at the weak sale is independent of its value at the GUT 
scale, provided the Yukawa couplings at the unification scale are large. For the fixed point 
or quasi-fixed point scenarios to be successful, it is necessary that these fixed points be 
stable!, I, §. 

Since supersymmetry necessciates the introduction of superpartners for all known 
particles in the SM (in addition to the introduction of two Higgs doublets), which trans- 
form in an indentical manner under the gauge group, we have additional Yukawa couplings 
in supersymmetric models which violate baryon number (B) or lepton number (L). In 
the MSSM a discrete symmetry called R-parity (Rp) is invoked to eliminate these B and L 
violating Yukawa couplings Pj. However, the assumption of Rp conservation at the level 
of MSSM appears to be ad hoc, since it is not required for the internal consistency of the 
model. Therefore, the study of MSSM, including R-parity violation, deserves a serious 
consideration. 

Recently attention has been focussed on the study of renormalization group evolution 
of Rp violating Yukawa couplings of the MSSM 0, and their quasi- fixed points. This 
has led to certain insights and constraints on the quasi-fixed point behavior of some of 
the Rp violating Yukawa couplings, involving higher generation indices. We recall that 
the usefulness of the fixed point and quasi-fixed point scenarios is the existence of stable 
infra-red fixed points. The purpose of this paper is to address the important question of 
the infra-red fixed points of supersymmetric models with Rp violation, and their stability. 
Our interest is in the structure of the infra-red stable fixed points, rather than the actual 
values of the fixed points. 

To this end we shall consider the supersymmetric standard model with the minimal 
particle content and with Rp violation, and refer to it as MSSM with R-parity violation. 
We begin by recalling some of the basic features of the model. The superpotential of the 
MSSM is given by 

W = flH,H2 + {K)abQiU'RH2 + {h,)abQlD'^H, + {hE)abLlE''^H,, (1) 



to which we add the L and B violating terms 

Wl = fi^L,H2 + ^XatcLlLlEl^ + X'^,^LlQlD% (2) 

Wb = iKtcDnDnUl, (3) 

respectively, as allowed by gauge invariance and supersymmetry. In Eq. (|1]), {hu)ab, (ho) at 
and {hE)ab are the Yukawa coupling matrices, with a, 6, c as the generation indices. The 
Yukawa couplings Xabc and A"j,^ are antisymmetric in their first two indices due to SU{2)i 
and SU{3)c group structure. Phenomenological studies of supersymmetric models of this 



type have placed constraints ||10l on the various couplings Xabc, X'^f^^ and A'^'^^, but there 
is still considerable room left. We note that the simultaneous presence of the terms in 
Eq. (0) and Eq. (^) is essentially ruled out by the stringent constraints ||Tl| implied by 



the lack of observation of nucleon decay. 

In addition to the dominant third generation Yukawa couplings ht = {hjj)^^-, hb = 
(^d)33 and hj. = (/i£;)33 in the super potential (|ip, there are 36 independent Rp violating 
couplings Xabc and A^^^ in Eq. (||), and 9 independent A"^^ in Eq. @. Thus, one would 
have to solve 39 coupled non-linear evolution equations in the L-violating case, and 12 
in the B-violating case, in order to study the evolution of the Yukawa couplings in the 
minimal model with Rp violation. In order to render the Yukawa coupling evolution 
tractable, we need to make certain plausible simplications. Motivated by the generational 
hierarchy of the conventional Higgs couplings, we shall assume that an analogous hierarchy 
amongst the different generations of Rp violating couplings exists. Thus we shall retain 
only the couplings A233, A333 and A233, and neglect the rest. We note that the Rp violating 
couplings to higher generations evolve more strongly because of larger Higgs couplings 
in their evolution equations, and hence could take larger values than the corresponding 
couplings to the lighter generations. Furthermore, the experimental upper limits are 
stronger for the Rp violating Yukawa couplings corresponding to the lighter generations. 

We shall first consider the evolution of Yukawa couplings arising from superpoten- 
tials ([^) and (^, which involve baryon number violation. The one-loop renormalization 
group equations for ht, hh, h^ and A233 (all others set to zero) are: 



16>r^5j^ = h {hi + 6hl + AJ + 2Ag3 - |gj - 3gl - ^gl 

2 dhr _ 7 foh2 , AU2 o„2 9„2 



For completeness we list the well-known evolution equations for the gauge couplings, 
which at one-loop order are identical to those in the MSSM, since the additional Yukawa 



coupling(s) do not play a role at this order: 

2 ddi _ 



IGtt 



Mf, , = 1,2,3, 



with hi = 33/5, 62 = 1, ^3 



(i(ln/i) 
3. With the definitions 



(5) 
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(6) 



and retaining only the SU{3)c gauge coupling constant, we can rewrite the renormaliza- 
tion group equations as {a^ = gl/i^lGn^)): 



dRf 



dt 



dRi) 
If 



asRt 

d^Rb 

dRr 

"df 



16 



6Rt - i?6 - 2R" 



16 

y 



b-A -Rt- QRb -Rr- 2i?" 

dsRr [h - ^Rb - ^Rr] , 



dR" 

"df 



daR" [(8 + 63) - 2Rt - 2Rb - QR"] 



(7) 

(8) 

(9) 

(10) 



where 63 = —3 is the beta function for g^ in the MSSM, and t = — In /i^. Ordering the 
ratios as Ri = {R", Rr,Rb, Rt), we rewrite the RG equations (|^ - (|1^) in the form ^ 



dRi 
~di 



asRi 



+ 63) - E ^^.^. 



'111 



where r^ = J2r 2 Cr, Cr is the QCD Casimir for the various fields {Cq = Cjj = Cjj = 4/3), 
the sum is over the representation of the three fields associated with the trilinear coupling 
that enters Ri, and S* is a matrix whose value is fully specified by the wavefunction 
anomalous dimensions. A fixed point is then reached when the right hand side of Eq. ([TT| ) 
is for all i. If we were to write the fixed-point solutions as R*, then there are two fixed 
point values for each coupling: R* = 0, or 



{n + h)-J2S^JR* 
j 

It follows that the non-trivial fixed point solution is 

R: = Y.iS-%ir, + b3] 



(12) 



(13) 



Since we shall consider the fixed points of the couphngs ht, hb and A233 only, we shall 
ignore the evolution equation @. However, the coupling h^- does enter the evolution 



Eq. (Q) of hb, but it can be related to hf, at the weak scale (which we take to be the 
top-quark mass), since 

hr{mt) = ^, 14 

rJT-v cos p 

and 

hr{mt) = ^, \ —hbirrit) = O.Qh{mt), (15) 

mb{mb) r]r 

where rjb gives the QCD or QED running [|T2| of the b-quark mass mb{fi) between ji = rrib 
and fj, = rrit (similarly for rj^-), and tan/5 = W2/'yi is the usual ratio of the Higgs vacuum 
expectation values in the MSSM, with v = {\/2Gf)~^^'^ = 246 GeV. The anomalous 
dimension matrix S can, then, be written as 



S=\ 2 6 + T] 1 , (16) 




where rj = h'^{int) / h\{mt) ^ 0.36 is the factor coming from Eq. (p!5|). We, therefore, get 
the following fixed point solution for the ratios: 

Rl ^ R"* = J!^±Z% ^ 0.76, 
^ 3(170 + 32r/) 

20 

fl. . r; = ^^±i!L ^ 0.12, 

^ * 170 + 32r7 

Since each of the -Rj's is positive, this is a theoretically acceptable fixed point solution. 

We next try to find a fixed point solution with R"* = 0, with Rb and Rt being given by 
their non-zero solutions. We need to consider only the lower right hand 2x2 sub-matrix 
of the matrix S in Eq. (|16D to obtain the fixed point solutions for Rb and Rt in this case. 
We then have 

Rl = R"* = 0, 
35 

^^-^-3(35Tm"°-''' ^''^ 

R* ^ R* = U^ + V) ^ Q 34^ 
^ * 3(35 + 6r]) 

This is also a theoretically acceptable solution, as all the fixed point values are non- 
negative. We must also consider the fixed point with Rl = 0, which is relevant for 
the low values of the parameter tan/3. In this case, we have to reorder the couplings as 
Ri = {Rb, R", Rt), so that we have the anomalous dimension matrix (in this case denoted 
as S) 

' 6 + r] 2 1 
^ = I 2 6 2 I . (19) 

1 2 6 



Since R^ = 0, we have to determine the non-zero fixed point values for R" and Rt only. 
For this we consider the lower right hand 2x2 submatrix of the matrix in (|T9|) to obtain 



iXn 



p* p* 

= R"* 



p* — D* 



19 ~ 79 

24 ~ ' 

- ^0.12. 



(20) 



which is an acceptable fixed point solution as well. Since there are more than one the- 
oretically acceptable IRSFPs in this case, it is important to determine which, if any, is 
more likely to be realized in nature. To this end, we must examine the stability of each 
of the fixed point solutions. 

The infra-red stability of a fixed point solution is determined by the sign of the eigen- 
values of the matrix A whose entries are {i not summed over) 



A 



u 






i]i 



f211 



where R* is the set of the fixed point solutions of the Yukawa couplings under consid- 
eration, and Sij is the matrix appearing in the corresponding RG equations (^TJ) for the 
ratios Ri. For stability, we require all the eigenvalues of the matrix Eq. (|21|) to have 
negative real parts (note that the QCD /?- function 63 is negative). Considering the fixed 
point solution ([T7|), the matrix A can be written as 



-4 


/ QRl 2RI 2RI \ 

2R* {6 + ri)R; R* , (22) 
\ 2i?3 ^3 QR*3 ) 


where R* are given in Eq. ([171). The eigenvalues of the matrix Eq. (^2]) are calculated to 

bp 


Ai = -1.6, A2 = -0.2, A3 = -0.2, (23) 


which shows that the fixed point ( 


TTD is an infra-red stable fixed point. We note that 



the eigenvalue Ai is larger in magnitude as compared to the other eigenvalues in ( pBj) 
indicating that the fixed point for A233 is more attractive, and hence more relevant. 



Next, we consider the stability of the fixed point solution (|TB|). Since in this case the 
fixed point of the coupling R"* = 0, we have to obtain the behaviour of this coupling 
around the origin. This behaviour is determined by the eigenvalue 



Ai 



1 



i=2 



p* 
Rj- 



fi + h 



(24) 



where ri = 2{C-fj+ Cjy) = 8, the Cs are the quadratic Casimirs of the fields occuring in 
the B- violating terms in the superpotential (^, and the S^j is the matrix (0), with the 
fixed points R*, i 



1, 2, 3 given by Eq. (18). Inserting these values in Eq. 



Ai 



385 + 76r/ 
9 (35 + 6ri) 



>0, 



we find 



(25) 



thereby indicating that the fixed point is unstable in the infra-red. The behaviour of the 
couphngs Rb and Rt around their respective fixed points is governed by the eigenvalues 
of the the 2x2 lower submatrix of the matrix A in Eq. ( p2D 



l/(6 + r/)/?* R; 

3 I Rl QRl 



(26) 



which we find to be 

A2 = -0.78, A3 = -0.56. (27) 

Although A2 and A3 are negative, because of the result (^), the fixed-point solution ([18|) 



is unstable in the infra-red. In other words, the Rp conserving fixed point solution ([18 
will never be achieved at low energies and must be rejected. 

Finally we come to the question of the stability of the fixed point solution (PP]). The 
behaviour of the coupling R^ around the origin is determined by the eigenvalue 



^'4 



i=2 



:28) 




where ri = 2(Cq + Cjy) = 16/3, and S is the matrix (|T^. Inserting these numbers, we 
find 

Ai = ^ ~ 0.21 > 0, (29) 

with the other two eigenvalues for determining the stability given by the eigenvalues of 
the matrix which is obtained from the lower 2x2 submatrix of the matrix 5* in (p^ . This 
submatrix can be written as 

(30) 

where R2 and -R3 are given by Eq. (|20D. The eigenvalues are 

A2 = -1.61, A3 = -0.22. (31) 

It follows, once again, that the fixed point solution given in ( pOD is not stable in the 
infra-red and is, therefore, never reached at low-energies. 

One may also consider the case where the couplings A233 and hb attain trivial fixed 
point values, whereas ht attains a non-trivial fixed point value. In this case we have R^ = 
Rl = 7/18, the well-known Pendleton-Ross |l|] top-quark fixed point of the MSSM. To 
study the stability of this solution in the present context, we must consider the eigenvalues 

O3 

where Si^ are read off from the matrix (|16D , which yields 

38 , _ 35 
-^^-27' ^^-54- 



Since the sign of each of Ai and A2 is positive, this solution is also unstable in the infra-red 
region. However, from our discussion of infra-red fixed point solution (0), it is clear that 
the Pendelton-Ross fixed point would be stable in case hh and A233 are small, though 
negligible at the GUT scale. In this case, these would, of course, evolve away from zero at 
the weak scale, though realistically they would still be small (but not zero) at the weak 
scale. Thus, the only true infra-red stable fixed point solution is the baryon number, and 
Rp, violating solution (p!?]). This is one of the main conclusions of this paper. We note 
that the value of R^ in {\[7t ) is lower than the corresponding value of 7/18 in MSSM with 
Rp conservation. 

It is appropriate to examine the implications of the value of ht{mt) predicted by our 
fixed point analysis for the top-quark mass. From (p!7D, and 0:3(7714) ~ 0.1, the fixed point 
value for the top- Yukawa coupling is predicted to be ht{mt) ~ 0.4. This translates into 
a top-quark (pole) mass of about nit — 70 sin /3 GeV, which is incompatible with the 
measured value [0 of top mass, rrit ~ 174 GeV, for any value of tan/3. It follows that 



the true fixed point obtained here provides only a qualitative understanding of the top 
quark mass in MSSM with Rp violation. 

We now turn to the study of the renormalization group evolution for the lepton number 
violating, and Rp, violating couplings in the superpotential (0). Here we shall consider 
the dimensionless couplings A233 and A333 only. The relevant one-loop renormalization 
group equations are: 

2 dhj- _ / 2 , .,2 , A\2 , Q\/2 



16vr^^^ = K [3ht + Ahl + 4A^33 + ?>\'i,^) , (32) 



dX') 

- '233 l^'V ^ ^^233 ^ ^^333 i ' 



167r^;^f^ = A233f4/7H4AL. + 3A'^ 
d[m jj,) 



levr^^ = A:,33 {hi + 6hl + hl + AL + 6A-33 - fgl) 



Defining the new ratios 



\2 \' 

p ^^233 p/ ^333 /QQ\ 

it — ^, rC — ^, lool 

93 m 



we may now rewrite the equations (p2]) as 



— = d^R [63 - 4i? - ?>R! - ARr] , (34) 



dR' 

asR' 



/16 



^^ — LV3 ' ^3) - i? - 6i?' - i?. - 6i?, - i?i , (35) 

dN 

" dsRr h -AR- 3R' - 3Rb - AR^] , (36) 



dt 



dRh 



dt 



tts-Rfe 



16 , 



a^Rt 



QR — Rt — QRb — Rt 

R — i?fe — QRt 

Ordering the ratios as Ri = {R, R', Rr, Rb, Rt), we can write these RG equations as: 

dRj 



dRt 

~~dt 



3 ^ 



dt 



a^Ri 



+ 63) - E ^^.^. 



(37) 
(38) 

(39) 



where r^ = Y1,r '^Cr-, with Cr denoting the quadartic Casimir of the each of the fields, 
the sum being over the representation of fields that enter i?j, and S fully specified by the 
respective wavefunction anomalous dimensions. It follows that there are two fixed point 
values for each coupling: i?* = 0, or the non-trivial fixed point solution 



R: = Y.is-%{r,+h,). 



(40) 



We shall be interested in the fixed-point solutions of the couplings A233, A333, hb, ht only, 
and shall not consider the h^- coupling. Therefore, we replace it, as we did earlier, by 
hr{mt) = O.Ghblmt) at the weak scale in the determination of the fixed point solutions 
(HOl). The anomalous dimensions matrix can then be written as: 



S 



1 


Vo 



6 + r/ 

6 + r/ 

1 



1 
1 
6/ 



(41) 



This leads to the fixed point values for the ratios: 



Rl 
ttr, = rt 






p* — p* 



R* = 0, 
315 + 19Ar] 
3661] — 315' 
140 



122r/ 
llOr]- 



-105 
105 



(42) 



366r7 — 315 

R'* < 0, and therefore, this fixed point solution is not an acceptable 
fixed point. We, thus, see that a simultaneous fixed point for the lepton number violating 
couplings A233, A333, and hb, ht does not exist. 

We now consider the two L- violating couplings separately, i.e., we shall take either 
333, or A333 -C A233, respectively. In the case when A333 is the dominant of the 



We note that R2 



A233 ^ A 



couplings, we order the couplings as Ri 
Eq. (EOl) for this case can be written as 



{R', Rb, Rt), so that the matrix S that enters 



6 6 + 7] 1 

S = \ 6 6 + r/ 1 

1 1 6 



(43) 



Since the determinant of this matrix vanishes, there are no fixed points in this case. We 
thus conclude that a simultaneous non-zero fixed point for the coupling A233, hb, ht does 
not exist. We note that the vanishing of the determinant corresponds to a solution with 
a fixed line or surface. 

If hh is small (e.g., for the case of small tan/?) we may reorder the couplings Ri = 
{Rb, R', Rt), and the matrix S, to find the fixed point solution 



p* — D* 
iXo = JX 



p* 



0, 

1 

3' 

r; = -. 
* 3 



(44) 



In order to study the stability of this solution, we must obtain the behaviour of the 
coupling Rl around the origin from the eigenvalue 



03 



"^SijR* - (ri + 63 
i=2 



where ri = 16/3. Inserting the relevant -R*'s into (^5|), we get 



(45) 



(46) 



from which we conclude that the fixed point Eq. (|^) will never be reached in the infra- 
red region. This fixed point is either a saddle point or an ultra-violet fixed point. We 
conclude that there are no non-trivial stable fixed points in the infra-red region for the 
lepton number violating coupling A333. 

Finally, we consider the case when A333 ^ A233. We find the fixed point solution 



Rl: 
p* 

Rt 



R* 



Rt 



R*t 



-315- 194r] 
12(35 + 6ri) ' 

35 

3(35 + 6r/)' 

7(5 + r/) 



(47) 



3(35 + 6r/)' 
which is unphysical. We, therefore, try a fixed point with Rl = 0. We find 



p* — p* 
p* — p* 



p* 



R*t 



= 0, 

-3/4, 
7/18, 



(4J 



which, again, is unphysical. We have also checked that: (1) trivial fixed points for A233 
and hb and the Pendleton- Ross type fixed point for the top-quark Yukawa coupling, or (2) 
trivial fixed points for A333 and hb and the Pendleton-Ross fixed point for the top-quark 
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Yukawa coupling, are, both, unstable in the infra-red region. We, therefore, conclude that 
there are no fixed point solutions for the lepton number violating coupling A233. 

To summarize, we have analyzed the one-loop renormalization group equations for 
the evolution of Yukawa couplings in MSSM with Rp violating couplings to the heaviest 
generation, taking into account B and L violating couplings one at a time. The analysis of 
the system with Rp, and the baryon number, violating coupling A233 yields the surprising 
and important result that only the simultaneous non-trivial fixed point for this coupling 
and the top-quark and b-quark Yukawa couplings ht and hi, is stable in the infra-red 
region. However, the fixed point value for the top-quark coupling here is lower than its 
corresponding value in the MSSM, and is incompatible with the measured value of the 
top-quark mass. The Rp conserving solution with A233 attaining its trivial fixed point, 
with ht and hi, attaining non-trivial fixed points, is infra-red unstable, as is the case for 
trivial fixed points for A233 and h^, with a non-trivial fixed point for ht. Our analysis 
shows that the usual Pendelton-Ross type infra-red fixed point of MSSM is unstable in 
the presence of Rp violation, though for small, but negligible, values of hb and A233 it 
could be stable. The system with L, and Rp, violating couplings does not possess a set 
of non-trivial fixed points that are infra-red stable. Our results are the first in placing 
strong theoretical constraints on the nature of Rp violating couplings from fixed-point 
and stability considerations: the fixed points that are unstable, or the fixed point that is 
a saddle point, cannot be realized in the infra-red region. The fixed points obtained in 
this work are the true fixed points, in contrast to the quasi-fixed points of [H, and serve 
as a lower bound on the relevant Rp violating Yukawa couplings. In particular, from our 
analysis of the simultaneous (stable) fixed point for the baryon number violating coupling 
A233 and the top and bottom Yukawa couplings, we infer a lower bound on A233 ~ 0.98. 



Note addded: After this paper was submitted for publication, another paper \Ti\ which 
considers the question of infra-red fixed points in the supersymmetric standard model 
with Rp violation has appeared. The fixed points for A333, A233 and ht, neglecting all 
other Yukawa couplings, is considered. Their results, where there is an overlapp, agree 
with ours. However, unlike in the present work, the stability of the fixed points has not 
been considered in [Q. 
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